BARI-TH/284-97 



The solar neutrino problem 
after three hundred days of data at SuperKamiokande 



G. L. Fogli, E. Lisi, and D. Montanino 

Dipartimento di Fisica and Sezione INFN di Bari, 
Via Amendola 173, 1-70126 Bari, Italy 



Abstract 

We present an updated analysis of the solar neutrino problem in terms 
of both Mikheyev-Smirnov-Wolfenstein (MSW) and vacuum neutrino os- 
cillations, with the inclusion of the preliminary data collected by the Su- 
perKamiokande experiment during 306.3 days of operation. In particular, the 
observed energy spectrum of the recoil electrons from neutrino scattering is 
discussed in detail and is used to constrain the mass-mixing parameter space. 
It is shown that: 1) the small mixing MSW solution is preferred over the large 
mixing one; 2) the vacuum oscillation solutions are strongly constrained by 
the energy spectrum measurement; and 3) the detection of a possible semi- 
annual modulation of the flux due to vacuum oscillations should require 
at least one more year of operation of SuperKamiokande. 

PACS number(s): 26.65.+t, 13.15.+g, 14.60.Pq 
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I. INTRODUCTION 



The solar neutrino problem namely, the deficit of the neutrino rates measured by 
the four pioneering solar neutrino experiments, Homestake p[, Kamiokande [Q], SAGE [Q, 
and GALLEX [§, as compared to the standard solar model predictions [^, represents one 
of the most convincing indications for new physics beyond the standard electroweak theory. 
On the one hand, the confidence in the performances of the four detectors has increased 
continuously, as a result of many careful experimental cross-checks and calibrations p|-p!0|. 
On the other hand, the reliability of the most refined solar evolution models (including 
light element diffusion) has been corroborated by the impressive agreement with a growing 
amount of helioseismological data [|ll|,|l^. Therefore, new neutrino physics (such as neutrino 
oscillations) appears to be a likely solution of the solar v deficit. 

However, the deficit of the neutrino event rate is a solar model dependent quantity. It 
is highly desirable, instead, to measure observables whose interpretation is not related to a 
prior knowledge of the absolute neutrino flux, such as shape deviations of energy spectra, or 
relative variations of the event rates during the time of the year. Such measurements have 
been pioneered by the Kamiokande experiment p], but with statistics too low to provide 
definitive indications. This experimental program is now being pursued with much higher 
statistics by the real-time, water-Cherenkov SuperKamiokande experiment ||13|| . 

The preliminary results after the first 306.3 days of operation of SuperKamiokande 
||I^^TB[ are already sufficiently accurate to provide interesting new insights into the solar 
neutrino problem, although further data are required to draw decisive conclusions. There- 
fore, we think it useful to present an updated analysis of the solar neutrino problem, includ- 
ing the most recent results from the four pioneering experiments [ 17,3,18, 19 1 and from the 
SuperKamiokande experiment after 306.3 days [p!^ -[T6|, that have become publicly available 
during the 1997 summer conferences. In particular, we discuss in detail the information 
provided by the energy spectrum of recoil electrons observed at SuperKamiokande, which is 
one of the most important solar model independent observables. The available experimental 
information is interpreted in the light of two-flavor neutrino oscillations, both in vacuum 
and in matter, according to the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism . 

The paper is organized as follows. In Sec. II we analyze the solar neutrino problem using 
the available solar model dependent information (namely, the absolute neutrino rates). In 
Sec. Ill we analyze the solar model independent information, and in particular the elec- 
tron energy spectrum measured by SuperKamiokande, in order to constrain the oscillation 
interpretation both in matter and in vacuum. We draw our conclusions in Sec. IV. Some 
technical aspects of the energy spectrum analysis are elucidated in the Appendix. 



II. SOLAR MODEL DEPENDENT INFORMATION 

In this section we report the most recent measurements of the solar neutrino rates, and 
compare them to the expectations of the 1995 Bahcall-Pinsonneault (BP95) standard solar 
model 0]. The observed deficit is interpreted in terms of MSW and vacuum two-family 
oscillations in the parameter space spanned by the neutrino squared mass difference bvn? 
and by the mixing angle Q. 
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A. The solar neutrino deficit 



Table I reports the neutrino event rates measured by the four pioneering solar neutrino 
experiments Homestake [|1^], Kamiokande P], SAGE ||18[, and GALLEX together with 



the neutrino flux measurement after 306.3 days at SuperKamiokande [|14|-[16[|. The ob 



served rates appear to be significantly smaller than the theoretical expectations of the BP95 
model. 

Figure 1 shows graphically the information reported in Table I, with the further inclu- 



sion of the correlations of theoretical uncertainties, which have been calculated as in p2 
In Fig. 1 the GALLEX and SAGE data have been combined (in quadrature) in a single 
(Gallium) result. The same has been done for the Kamiokande and SuperKamiokande data 
(K+SuperK). In the combination, asymmetric errors have been conservatively symmetrized 
to the largest one. The 99% C.L. experimental and theoretical ellipses appear to be distinctly 
separated in all the planes charted by any two experiments. Notice that the (strongly corre- 
lated) theoretical errors are rather large as compared to the experimental errors. Therefore, 
as far as the total neutrino rates are concerned, further reductions of the experimental un- 
certainties are not expected to change significantly the current picture of the solar v deficit, 
while improvements in the theoretical predictions would have a greater impact. 



B. MSW and vacuum oscillation fits 

As is well known, a viable explanation of the solar neutrino deficit is represented by 
matter-enhanced (MSW) neutrino oscillations (see, e.g., p3|-p5|) or, alternatively, by vacuum 
neutrino oscillations (see, e.g., p6| , p5| ). Assuming for simplicity oscillations between two 
families, the analysis involves only two parameters, the neutrino squared mass difference 
bvn? and the mixing angle Q. 

Figure 2 shows the results of our x^-analysis of the data reported in Table I, assuming 
MSW oscillations. Notice that this fit includes only the deficit data, and not the solar 
model independent information provided by the energy spectra or the night-day asymmetry 
(that will be considered separately in Sec. III). The usual small and large mixing angle 
solutions appear to be slightly more constrained than in previous fits (see, e.g., p^), mainly 
as a result of the new SuperKamiokande data. The absolute minimum of (Xmin = 1-0) 
is reached at (5m^, sin^2^) = (9.6 x 10^^ eV^, 4.8 x 10"'^); the secondary minimum (xLc = 
2.45) is reached at (^m^, sin^ 2Q) = (L5 x IQ-^ eV^ 0.58). The 90, 95, and 99% C.L. allowed 
regions coincide with the contours at ~ Xmin = 4.61, 5.99, and 9.21, respectively. 

Figure 3 shows the analogous results for the vacuum oscillation analysis. We find the 
minimum value xLin = 4.0 at (5m^, sin^ 29) = (8.2 x 10~^^ eV , 0.87). A comparison of Xmm 
seems to indicate a preference of the MSW oscillation scenario over the vacuum oscillation 
one. However, it should be noted that, in vacuum oscillation fits, the value of Xmin is very 
sensitive to small shifts in the central values of the experimental results and thus it might 
change significantly with new data (while it turns out to be more stable in MSW fits). 

A few technical remarks are in order. In our oscillation analyses, the GALLEX and SAGE 
data have been combined in quadrature. However, the Kamiokande and SuperKamiokande 
data have been fitted separately, since these two experiments have different thresholds and 
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energy resolutions. In particular, according to [16|, we use for SuperKamiokande a threshold 



of 6.5 MeV for the total electron energy E^, and a Gaussian energy resolution function 
with a la width of ±15% at T = 10 MeV {T = — rrif,), scaling as 1/VT at different 
electron kinetic energies. The input neutrino fluxes have been taken from 0, and the 
corresponding theoretical uncertainties have been included in the covariance matrix as 
in [^. Asymmetric errors have been conservatively symmetrized to the largest one. The 
®B neutrino energy spectrum has been taken from [^.Q Concerning the neutrino cross 
sections, we use the updated calculations of for the chlorine experiment, and of ||2^ 



for the water-Cherenkov experiments. In the MSW calculations, the neutrino production 
regions and the electron density profiles have been taken from p. The Earth regeneration 
effect is included analytically as in ||29| for each detector latitude. In the vacuum oscillation 
calculations, the time average over the year is performed through the analytical approach 
of [Q. In conclusion. Figs. 2 and 3 represent state-of-the-art results in the field of solar 



neutrino oscillations. 



III. SOLAR MODEL INDEPENDENT INFORMATION 

In this section we analyze the solar model independent information about the electron 
energy spectrum and the time variation of the neutrino rate measured by SuperKamiokande. 
Since these data are already more accurate than the corresponding Kamiokande ones, we 
do not include the latter in our analysis. 

We recall that the SuperKamiokande experiment has collected a total of 4395 solar neu- 
trino events above threshold {E^ > 6.5 MeV) during 306.3 days of operation with a fiducial 
volume of 22.5 kton ||I^-[16|, corresponding to an observed rate of 0.6377 events/day/kton. 



Given the quoted deficit factor of 0.3685 Ijl^ (see also Table I), the expected rate (BP95 
model, no oscillation) is equal to 1.730 events/day/kton. 



A. Energy spectrum 



The SuperKamiokande experiment has measured the energy spectrum of recoil electrons 
from u-e scattering. The spectrum is given [|l4|-0 in the form of a 16-bin histogram, 
each bin presenting the ratio between the experimental rate and the theoretical rate. Since 
the theoretical spectrum has not been explicitly presented, we rely upon our calculations, 
using for SuperKamiokande the same inputs (^B u spectrum, cross section, energy resolution 
and threshold) described in the previous section, and normalizing the results to the total 
expected rate of 1.730 events/day/kton. The global information about the energy spectrum 
is reported in Table II and in Fig. 4. 



The B spectrum shape errors evaluated in |27| have not been included in the fits of Figs. 2 and 
3, since their effect is much smaller than the uncertainty in the absolute neutrino flux. However, 
such shape errors will be included in the analysis of the energy spectrum performed in Sec. Ill and 
in the Appendix. 
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Table II shows the relevant spectral data in each of the 16 bins (numbered in the first 
column). As shown in the second column, all bins have a 0.5 MeV width, with the exception 
of the last one, which collects all events between 14 and 20 MeV (total energy). The third 
and fourth column report the average total energy (in MeV) and the expected event rate 
rii (in events/day/kton/bin) for each bin, in the absence of oscillations (our calculation). The 
sum of the entries in the fourth column gives the expected rate of 1.730 events/day/kton. 
The fifth column reports the ratio between the experimental and theoretical rate in each 
bin. The upper and lower errors of rj are given in the following columns. The sixth and 
seventh columns report the quadratic sum of statistical and uncorrelated systematic errors 
(cr-). Finally, the correlated systematic errors {a") are given in the last two columns. The 
data in columns 6-9 have been graphically reduced from the plots shown in 0-|l6| and thus 



may be subject to slight inaccuracies. The numbers are given with three decimal places just 
to avoid further errors due to truncation and round-off. The asymmetry between upper and 
lower errors is mainly due to two slightly different energy calibration procedures |]TB| . Shifts 



in the absolute energy scale represent the most important source of systematic uncertainties, 



as was emphasized earlier in 1 3 1,32 



Figure 4 shows the spectral information in a different way. The solid curve represents 
our calculation of the electron spectrum, using the best-fit neutrino spectrum given in 
27|. The shape of the neutrino spectrum is subject to some uncertainties, that have 



been carefully evaluated in |2^. Therefore, we have also used the ±3(T deviated neutrino 
spectra reported in ||2^ to calculate the effect on the electron spectrum (dotted curves). 
Notice that in Fig. 4 both the central and the ±3a deviated theoretical electron spectra are 
renormalized to the same total area as the experimental spectrum (0.6377 events/day/kton). 
The black circles represent the experimental spectrum, with horizontal error bars spanning 
the bin widths, and vertical error bars representing the la total experimental uncertainties 
(quadratic sum of statistical and systematic errors). It can be seen that the slope of the 
experimental spectrum is slightly more gentle than the theoretical spectrum, the rate ob- 
served at low (high) energies being somewhat suppressed (enhanced) with respect to the 
expectations. This indication, if confirmed with significantly smaller errors, would repre- 
sent unmistakable evidence for a new, energy- dependent process (such as flavor oscillations) 
affecting solar neutrinos along their path to the Earth. 

Since the observed spectrum deviates only slightly from the expected one, it is reasonable 
to try to "summarize" the spectral information in a single parameter, related in some way to 



the slope deviation. We adopt the approach advocated in Refs. [^,^,^ , where the spectral 
deformations were expanded in a series of deviations of the spectral moments (the mean, the 
variance, etc.) from their standard values. Given the present experimental uncertainties, 
it is sufficient to study the deviations of the first moment, namely, of the average kinetic 
energy (T) of the electrons, from its standard (no oscillation) value {T)q. This reduction 
procedure has the practical advantage that a single parameter (A(T)) is used in the fits, 
instead of 16 bins with correlated errors. 

Since the value of (T) has not been reported in IT^HTH], we have estimated it from the 



energy spectrum information presented in Table II and in Fig. 4. The reader is referred to 
the Appendix for a detailed derivation. Our result for the fractional deviation (T)/ (T)o — 1 
is: 
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X 100 = 0.991-^^ , (1) 

where the errors represent the total uncertainties at la. This result represents the starting 
point of our analysis of the energy spectra deviations due to oscillations. It can be noticed 
that the observed value of (T) is about la higher than the expectations, with an upper error 
larger than the lower one. Therefore, negative deviations of (T) (i.e., negative "tilts" of the 
spectrum slope) will be much more constrained than positive deviations. 



B. Time variations 

Possible variations of the neutrino rate during the night would represent evidence for the 
a regeneration of z/g's in the Earth matter, as expected within the MSW scenario (see, e.g., 
29| and references therein). The SuperKamiokande experiment has found no such effect 



within the present uncertainties. In fact, the measured asymmetry between nighttime (N) 
and daytime (D) neutrino rates is |I^-p!^ 



N - D 

^^j^ = 0.017 ± 0.026(stat.) ± 0.017(syst.) . (2) 

No evidence for variations of the neutrino rate during the time of the year has been found 
so far Semiannual modulations of the signal (in addition to the trivial 1/L^ geometrical 
variations) would indicate the presence of neutrino oscillations in vacuum (see, e.g., |]30|j33[| 
and references therein). A Fourier analysis of the signal expected in SuperKamiokande 
shows that this experiment is sensitive only to the first harmonic /i which, in the absence 
of oscillations, is equal in value to the Earth's orbit eccentricity e = 0.0167 (see for 
details). The purely statistical error af of the difference /i — e is given by [^] 



where Ns and Nb are the number of signal and background events, respectively. 

The numbers Ns and Nb depend on the cut applied to the fiducial angle in the direction 
of the sun (6'sun)- Using Eq. (H) and the spectrum of Ns + A''^ in terms of cos6'sun as 
reported in [p!^ -[l6|, we derive the values of aj shown in Table III. The dependence of aj on 
the cuts applied to cos^sun turns out to be weak, since the effect of a lower signal tends to 
be compensated by a a better signal-to-background ratio, and vice versa. On the basis of 
the results reported in Table III, we assume 

af - ±0.020 (1 yr) (4) 

as the typical statistical error of the first Fourier harmonic, after one year of operation at Su- 
perKamiokande. The errors corresponding to 2 and 4 years of operation are simply obtained 
with rescaling factors of 1/a/2 and 1/2, respectively (assuming the current background event 
rate). Possible systematic errors that might affect the determination of /i are not considered 
here. 
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C. Implications for oscillations in matter 



Figure 5 illustrates the implications of the electron energy spectrum measured by Su- 
perKamiokande [as summarized by the datum of Eq. (|l|)] for the MSW scenario. The upper 
panel shows curves at constant values of the fractional deviation of (T) expected in the 
presence of MSW oscillations. Similar curves were discussed earlier in p9[ (with a different 
energy threshold) and in (without the Earth regeneration effect). Superposed are the 
small and large mixing angle solutions found in Fig. 2 (at 95% C.L.). The central value 
in Eq. (|l|) is intriguingly close to the values of A{T) / {T)q within the small mixing angle 
solution, which thus appears slightly favored as compared to the large mixing one. Unfortu- 
nately, the errors in Eq. (|l]) are still large, so that only excluded regions can be meaningfully 
derived at present. The lower panel shows the regions excluded at 2, 3, and 4 standard 
deviations by Eq. (|l|). For instance, in the regions excluded at the 2a level, the theoretical 
value of A(T)/(T)o is either smaller than —0.93% or greater than +6.05%. The excluded 
regions are similar to those shown by the SuperKamiokande Collaboration |T^-|T6[. 

Figure 6 illustrates the implications of the night-day asymmetry datum of Eq. |^ for the 
MSW scenario. The upper panel shows iso-lines of the night-day asymmetry (percentage) 
expected in the presence of oscillations. The lower panel shows the regions excluded at 
2, 3 and 4 standard deviations by the datum of Eq. |^. The lower part of the large angle 
solution, where the asymmetry is predicted to be ~ 10%, appears to be disfavored by the 
experimental results. 

We do not attempt here a global x^-combination of all the SuperKamiokande data (total 
rate, energy spectrum, and night-day variations), since these pieces of information are not 
independent; in fact, they are just different "projections" of the double differential spectrum 
of events as a function of time and energy. A proper analysis of the whole SuperKamiokande 
data will be possible when such spectrum and its uncertainties will become available. 

In conclusion, the present data seem to favor the small angle solution with respect to the 
large mixing one, because: 1) The global value of is lower at the small mixing solution; 
2) The energy spectrum shows a slight deviation consistent with small mixing; 3) Part of 
the large mixing angle solution is disfavored by the nonobservation of a large night-day 
asymmetry. 



D. Implications for oscillations in vacuum 

Figure 7 illustrates the implications the energy spectrum measurement [as summarized 
by the datum in Eq. (|^)] for the vacuum oscillation scenario. The upper panel shows curves 
at constant values of the fractional deviation of (T) expected in the presence of vacuum 
oscillations. Similar curves were discussed earlier in |Q (with a different energy threshold for 
SuperKamiokande). Superposed are the multiple solutions found in the vacuum oscillation 
fit of Fig. 3 (at 95% C.L.). The solutions at "low ^m^" {5m'^ ~ 6-7.5 x lO'^^ eV^) predict 
positive deviations of (T), in agreement with the datum of Eq. ([^). The "high 5m^" solutions 
predict negative deviations, and thus are strongly disfavored, as shown in the lower panel of 



Fig. 7. Similar results have been discussed by the SuperKamiokande Collaboration |[l^-[16 
and are consistent with earlier Kamiokande spectrum fits |2 
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Figure 8 illustrates the sensitivity of the SuperKamiokande experiment to semiannual 
modulations of the neutrino rate induced by vacuum oscillations. The solid curve represents 
the deviations of the first Fourier coefficient /i from its standard value e (see |]3^), plotted as 
function of 5rn? in the most favorable condition (maximal oscillation amplitude, sin^ 26 = 1). 
The largest vertical error bar corresponds to the ±lcr statistical error estimated in Eq. (^ 
for one year of operation of SuperKamiokande. It can be seen that the 1 yr error bars are as 
large as the range spanned by variations of /i — e, which are thus practically undetectable at 
present. However, the error bars after 2 and 4 years of data taking should be reduced enough 
to allow some sensitivity to the largest predicted variations. Of course, improvements in the 
energy threshold or in the background rate might reduce the uncertainties more rapidly than 
the time sequence shown in Fig. 8. 

By comparing Figs. 7 and 8, it can be noticed that the deviations of (T) and of /i have 
the same sign (both positive or negative) at equal values of 5im?. Such correlation between 
energy spectrum deviations and time variations of the signal can be traced to the L/E 



dependence of the vacuum oscillation probability, as recently emphasized in [34 



In conclusion, the SuperKamiokande energy spectrum data constrain rather strongly the 
vacuum oscillation solutions, allowing only the restricted range 5m? ~ 6-7.5 x 10~^^ eV^ (at 
large mixing). The present experimental sensitivity does not allow detection of semiannual 
modulations of the neutrino rate due to oscillations. Detection at greater than la level 
should require at least one more year of operation, unless the uncertainties get reduced 
further by improving the energy threshold or the background rate, relative to the present 
levels. 



IV. SUMMARY AND CONCLUSIONS 

We have performed a thorough analysis of the solar neutrino problem in two scenarios: 
MSW and vacuum oscillations (between two neutrino families). We have used the most 
recent available data, including the preliminary results of the SuperKamiokande experiment 
after 306.3 days of operation. Particular attention has been devoted to the treatement of the 
new, solar-model independent data, namely, the energy spectrum information and the time 
variations of the signal in SuperKamiokande. It has been shown that: 1) the small mixing 
angle solution is preferred over the large angle one; 2) the vacuum oscillation solutions 
are strongly constrained by the energy spectrum measurements; and 3) the detection of 
semiannual modulations of the neutrino flux in SuperKamiokande should require at least 
one more year of operation. 
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APPENDIX: EVALUATION OF (T) FROM THE SUPERKAMIOKANDE 

SPECTRUM 



In this Appendix we give the detailed derivation of the result reported in Eq. (|l|). We 
make use of the preliminary SuperKamiokande spectrum data after 306.3 days of operation 



l^-p]6[, as reported in Table II and shown in Fig. 4. 

The average value of the total electron energy E^. from the unbinned theoretical spectrum 
n{Eg) (solid curve in Fig. 4) is given by 

f.20MeV 



dEe n{Ee) Ee 

(Ee)theo (unbinned) = ^^^^y = 8.536 MeV . (Al) 

dEe n{Ee) 



/6.5 MeV 

Using instead the binned theoretical spectrum rii from Table II one obtains: 

y^l6 77'. 

(Ee)theo(binned) = ' ' = 8-525 MeV . (A2) 

The difference between the values in Eqs. ( [Al| ) and ( |A2| ) is only ~ 0.1%, indicating that 
the effect of binning is not important in the evaluation of the average energy. We will 
anyway attach to (-E'e)theo a ±0.010 MeV "binning error" that, as we shall see, is an order 
of magnitude smaller than the present experimental uncertainties. Another small error is 
induced by the neutrino spectrum uncertainty (as shown by the dotted curves in Fig. 4), 
that we evaluate as A(i?e)theo = ±0.016 MeV at the la level. 

The average value of Ee derived from the (binned) experimental spectrum is given by 

(^e)exp = ffe ' ' ' = 8-604 MeV , (A3) 

where r, and Ei are given in Table II. In general, the squared error a^^^ associated to (-Ee)exp 
reads 

^{E) = 1^2^ ]r P^J ^3 , (A4) 

i=l j=l ^' I ^' j 

where cxj represents the error of and pij is the correlation matrix. The partial derivatives 
are easily obtained from Eq. ( |A3| ) : 

The correlation coefficients for the statistical and uncorrelated systematic errors are given 
by Pij = Sij, so that from Eq. (|A4|) one obtains 



o"(£;)(stat. ± uncorr. syst.) = ± 



16 

E 

1=1 



exp / 



<9r,; 
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where the erf s are taken from Table II, upper and lower errors being propagated separately. 
The remaining systematic errors are fully correlated [Q, i.e., pij = 1, so that 



a,,, (corr. syst.) = ± g a'l = | +°- (MeV) , (A7) 

where the erf's are given in the last two columns of Table II. 

Taking the difference between Eqs. and and including all the uncertainties, 

one obtains 

(^)e.p - (^)theo = 0.079 i[J:Ji 1°:^^° ± o.oio ± o.oie (MeV) , (as) 

where the first error is due to statistics and to uncorrelated systematics, the second to 
correlated systematics, the third to the electron spectrum binning, and the fourth to the 
neutrino spectrum uncertainties. In conclusion, passing from total to kinetic energies, it 
follows that 



(T) 



theo 



0.991^1^ % , (A9) 



having added in quadrature the la independent errors in Eq. ( [ASI ). This is the final result 
anticipated in Eq. (|l|). 
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TABLES 



TABLE L Neutrino event rates measured by solar neutrino experiments, and corresponding 
predictions from the BP95 standard solar model. The quoted errors are at la. 

Experiment Ref. Data Theory Units 

Homestake Q 2.54 ± 0.16 ± 0.14 9.31};^ SNU 

Kamiokande [|] 2.80 ± 0.19 ± 0.23 ^■^'^-ill 10^ cm'^ s"^ 

SAGE H 73±8.5tli 137t? SNU 

GALLEX H 76.2 ±6.5 ±5 137t? SNU 

SuperKamiokande Q 2.44 ± 0.06l[j:^^ 6.62^?;?^ 10^ cm-^ g-^ 



TABLE II. The electron energy spectrum at SuperKamiokande. First three columns: sequen- 
tial number, energy range, and average energy for each bin. Fourth column: expected number of 
events rii per bin per day per kiloton without oscillations (our calculation). Fifth column: ratio of 
measured to expected rates in each bin. Sixth and seventh columns: Upper and lower la errors 
of Tj from statistical and uncorrelated systematic uncertainties. Eighth and ninth columns: Upper 
and lower la errors of from correlated systematic uncertainties. The numbers in columns 6-9 
have been graphically reduced from the plots shown in [p!4|-p!6|. 



i 


Range ^ 


Ei ^ 


ni ^ 




±0"- (stat.+unc.) 


±(T-' (corr. 


syst.) 


1 


[6.5, 7] 


6.747 


0.309 


0.346 


+0.046 


-0.032 


+0.010 


-0.004 


2 


[7, 7.5] 


7.244 


0.271 


0.341 


+0.045 


-0.032 


+0.013 


-0.006 


3 


[7.5, 8] 


7.742 


0.232 


0.350 


+0.043 


-0.036 


+0.016 


-0.008 


4 


[8, 8.5] 


8.239 


0.198 


0.400 


+0.062 


-0.037 


+0.025 


-0.010 


5 


[8.5, 9] 


8.742 


0.165 


0.327 


+0.064 


-0.039 


+0.026 


-0.012 


6 


[9, 9.5] 


9.240 


0.135 


0.373 


+0.075 


-0.046 


+0.040 


-0.014 


7 


[9.5, 10] 


9.736 


0.109 


0.382 


+0.092 


-0.049 


+0.050 


-0.016 


8 


[10, 10.5] 


10.239 


0.087 


0.373 


+0.101 


-0.053 


+0.058 


-0.017 


9 


[10.5, 11] 


10.737 


0.066 


0.356 


+0.115 


-0.059 


+0.069 


-0.023 


10 


[11, 11.5] 


11.234 


0.050 


0.372 


+0.138 


-0.062 


+0.089 


-0.026 


11 


[11.5, 12] 


11.736 


0.036 


0.392 


+0.170 


-0.072 


+0.110 


-0.036 


12 


[12, 12.5] 


12.234 


0.026 


0.402 


+0.209 


-0.088 


+0.137 


-0.042 


13 


[12.5, 13] 


12.730 


0.017 


0.359 


+0.229 


-0.098 


+0.141 


-0.046 


14 


[13, 13.5] 


13.232 


0.012 


0.422 


+0.238 


-0.123 


+0.192 


-0.063 


15 


[13.5, 14] 


13.730 


0.007 


0.627 


+0.506 


-0.183 


+0.337 


-0.101 


16 


[14, 20] 


14.817 


0.010 


0.493 


+0.510 


-0.170 


+0.378 


-0.114 



'^Units: MeV (total electron energy). 
^Units: events/day/kton/bin. 
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TABLE III. Number of signal events Ns and signal-to-background ratio Ns/Nb for various 
cos ^sun thresholds (as derived from the cos ^sun 

event spectrum presented in |p^-|l^), together 
with the corresponding uncertainty crj of the deviation of the first Fourier coefficient /i from its 
standard value e [see Eq. dH)]. 

cos^sun Ns Ns/Nb 

> 0.5 ~ 4100 ~ 0.30 0.023 

> 0.6 ~ 3870 ~ 0.35 0.022 

> 0.7 ~ 3620 ~ 0.38 0.021 

> 0.8 ~ 3270 ~ 0.59 0.020 

> 0.9 - 2530 ~ 0.92 0.020 
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FIGURES 



FIG. 1. The current state of the solar neutrino deficit. The ellipses represent the regions 
allowed at 99% C.L. by the present solar neutrino data (dashed lines) and by the standard solar 
model (SSM) of Bahcall-Pinsonneault [^. The coordinates are the chlorine (CI), gallium (Ga), and 
water-Cherenkov [normalized (K+SuperK)/SSM] neutrino rates. 



FIG. 2. The MSW small and large mixing angle solutions to the solar neutrino problem, as 
obtained by a fit to the total rates only (energy spectrum and night-day data not included). 



FIG. 3. The vacuum oscillation solutions to the solar neutrino problem, as obtained by a fit 
to the total rates only (energy spectrum and night-day data not included). The range of the 
(logarithmic) vertical scale is (0.5-2) x 10^^'^ eV^. 



FIG. 4. The electron recoil energy spectrum at SuperKamiokande. The dots represent the 
measured rates in each energy bin, with vertical la error bars given by the quadratic sum of 
statistical and systematic errors. The horizontal bars span the bin widths. Also shown is the 
standard spectrum expected in the absence of oscillations (solid line), together with the spectra 
obtained by assuming ib3cr deviations of the standard ®B neutrino spectrum [27|. The data are 
taken from The theoretical spectra refer to our calculations, renormalized to give the same 

area as the experimental histogram. 



FIG. 5. Constraints on the MSW solutions as derived by the measured value of the average 
electron kinetic energy (T). Upper panel: Theoretical expectations for the fractional shift of (T) 
from its standard (no oscillation) value {T)o- Lower panel: regions excluded at 2, 3, and 4 standard 
deviations by the SuperKamiokande determination of (T). 



FIG. 6. Constraints on the MSW solutions as derived by the measured value of the asymme- 
try of nighttime (-Rn) and daytime (-Rd) event rates. Upper panel: theoretical expectations for 
(-Rn — Rd)/{Rn + Rd)- Lower panel: Regions excluded at 2, 3, and 4 standard deviations by the 
SuperKamiokande data. 



FIG. 7. Constraints on the vacuum oscillation solutions as derived by the measured value of 
the average electron kinetic energy (T). Upper panel: Theoretical expectations for the fractional 
shift of (T) from its standard (no oscillation) value (r)o. Lower panel: Regions excluded at 2, 3, 
and 4 standard deviations by the SuperKamiokande determination of (T). 



FIG. 8. Fourier analysis of time variations induced by vacuum oscillations. Solid curve: ex- 
pected deviations of the first Fourier coefficient /i from its standard value (equal to Earth orbit ec- 
centricity e) for sin^ 20 = 1. Error bars: Our estimate of the statistical errors of SuperKamiokande 
(including background) after 1, 2, and 4 years of data taking. See the text and |33| for details. 
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FIG. 1. The current state of the solar neutrino deficit. The ehipses represent the regions 
ahowed at 99% C.L. by the present solar neutrino data (dashed lines) and by the standard solar 
model (SSM) of Bahcall-Pinsonneault The coordinates are the chlorine (CI), gallium (Ga), and 
water-Cherenkov [normalized (K+SuperK)/SSM] neutrino rates. 
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MSW solutions 

Chlorine + Gallium + Kamloka + SuperKamioka (306.3 d) total rates 
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FIG. 2. The MSW small and large mixing angle solutions to the solar neutrino problem, as 
obtained by a fit to the total rates only (energy spectrum and night-day data not included). 
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FIG. 3. The vacuum oscillation solutions to the solar neutrino problem, as obtained by a fit 
to the total rates only (energy spectrum and night-day data not included). The range of the 
(logarithmic) vertical scale is (0.5-2) x 10~^° eV^. 
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FIG. 4. The electron recoil energy spectrum at SuperKamiokande. The dots represent the 
measured rates in each energy bin, with vertical la error bars given by the quadratic sum of 
statistical and systematic errors. The horizontal bars span the bin widths. Also shown is the 
standard spectrum expected in the absence of oscillations (solid line), together with the spectra 
obtained by assuming ztSu deviations of the standard neutrino spectrum [27|. The data are 
taken from p!^-]l^. The theoretical spectra refer to our calculations, renormalized to give the same 
area as the experimental histogram. 
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SuperKamiokande, energy spectrum 
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FIG. 5. Constraints on the MSW solutions as derived by the measured value of the average 
electron kinetic energy (T). Upper panel: Theoretical expectations for the fractional shift of (T) 
from its standard (no oscillation) value {T)q. Lower panel: regions excluded at 2, 3, and 4 standard 
deviations by the SuperKamiokande determination of (T). 
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FIG. 6. Constraints on the MSW solutions as derived by the measured value of the asymme- 
try of nighttime (-Rn) and daytime (Rb) event rates. Upper panel: theoretical expectations for 
(i?N — ^d)/(^n + Rd)- Lower panel: Regions excluded at 2, 3, and 4 standard deviations by the 
SuperKamiokande data. 
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SuperKamiokande, energy spectrum 




FIG. 7. Constraints on the vacuum oscillation solutions as derived by the measured value of 

the average electron kinetic energy (T). Upper panel: Theoretical expectations for the fractional 
shift of (T) from its standard (no oscillation) value (T)o. Lower panel: Regions excluded at 2, 3, 
and 4 standard deviations by the SuperKamiokande determination of (T). 
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SuperKamiokande, vacuum oscillations 
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FIG. 8. Fourier analysis of time variations induced by vacuum oscillations. Solid curve: expected 
deviations of the first Fourier coefficient /i from its standard value (equal to Earth orbit eccentricity 
e) for sin^ 26 = 1. Error bars: Our estimate of the statistical errors of SuperKamiokande (including 
background) after 1, 2, and 4 years of data taking. See the text and [^] for details. 
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